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Passive  ranging  with  the  SNAP  model: 
A  performance  study 

E.J.  Sullivan  and  K.  Rameau 


Abstract:  A  passive  tanging  technique  which  uses  SNAP,  the  SACLANT 
Centre's  normal-mode  model,  is  studied.  An  examination  is  made  of  the 
sensitivity  of  the  technique  to  errors  in  the  assumed  sound  velocity  profile 
and  tilting  of  a  ISA  vertical  receiving  array.  Also,  the  effect  of  both  white 
and  coloured  noise  is  investigated.  Results  indicate  that  for  the  scenario 
considered  (summer  profile  in  the  Mediterranean  with  a  water  depth  of 
130  m  and  source  ranges  of  10  and  25  km),  the  method  can  be  quite  robust 
to  sound  velocity  profile  errors  or  the  order  of  several  m/s  and  array  tilts 
on  the  order  of  1.5*.  Solutions  were  obtained  with  input  signal-to-noise 
ratios  as  low  as  - 13  dB. 

Keywords:  array  tilt  o  bottom  properties  o  coloured  noise  o  depth 
o  gaussian  noise  o  inverse  problem  o  maximum  likelihood  □  passive 
ranging  o  propagation  o  sensitivity  o  signal-to-noise  ratio  o  SNAP  o 
sound  velocity  profile  o  white  noise  . 
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1.  Introduction 

Given  a  sufficiently  accurate  oceanic  propagation  model  and  a  sufficiently  accurate 
set  of  parameters  for  this  model,  one  can  solve  either  of  the  following  two  problems: 
(1)  given  enough  information  concerning  one  or  more  acoustic  sources,  find  the 
resulting  acoustic  field  anywhere  in  the  ocean;  (2)  given  a  sufficient  sample  of 
the  acoustic  field  in  the  ocean,  find  the  position  of  the  source(s).  Problem  (1) 
is  commonly  referred  to  as  the  forward  problem  whereas  problem  (2)  has  become 
known  as  the  inverse  problem.  In  this  report  we  are  concerned  with  the  inverse 
problem  and  more  specifically  the  effects  of  errors  in  the  input  parameters  on  such 
a  problem. 


2.  Theory 


Our  particular  inverse  problem  is  defined  as  follows.  Given  the  measurements  from 
a  set  of  equally- spaced  hydrophones  on  a  vertical  array,  find  the  range  to  the  narrow- 
band  acoustic  point  source.  The  approach  is  based  on  the  fact  that  the  normal-mode 
model  of  propagation  permits  a  set  of  linear  equations  to  be  written  that  can  be 
inverted.  For  fixed  frequency,  range  and  water  density,  the  normal-mode  model 
expresses  the  far-field  acoustic  pressure  at  range  r  and  depth  I  as 


f’(r.r) 


V”  #n(*?#n(^o)  -q.r+it.r 

2-  ji r 

n=l  V*" 


(1) 


Here,  <t>n(z)  is  the  nth  modal  function  at  depth  z,  z0  is  the  depth  of  the  source,  an 
is  the  loss  function  for  the  nth  mode,  k„  is  the  horizontal  wave  number  for  the  nth 
mode,  and  /V  is  the  total  number  of  modes 

A  more  comprehensive  discussion  of  Eq.  (1)  can  be  found  in  [1,  2].  For  our  pur¬ 
poses  it  is  sufficient  to  say  that  Eq.  (1)  represents  a  range-independent  model  and 
therefore  the  modal  function  loss  factor  and  the  horizontal  wave  number  all  can 
be  computed  with  knowledge  of  the  sound  velocity  profile  (SVP),  water  depth  and 
bottom  properties  but  without  knowledge  of  the  range.  Furthermore,  since  the 
modal  functions  are  real,  the  range  information  is  contained  in  the  phase  only,  and 
therefore  the  source  depth  need  not  be  known. 

Defining 


Mm,n  —  $n{  zm  )/  \/km, 

A„  =*„<z,)e“»r, 

=  A-e*-', 


-  I  - 


(2a) 

(2b) 

(2c) 
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Eq.  (1)  can  be  written  as 

Pm  —  Ainj.nXn-  (3) 

Inversion  of  Eq.  (3)  yields  the  N  functions  which,  in  principle,  can  be  solved  for 
the  range.  With  the  addition  of  noise,  Eq.  (3)  becomes 

Pm  =  +  cm,  1  <  n  <  N, 

1  <  m  <  Af,  (4) 

where  A/  is  the  number  of  hydrophones  and,  as  defined  previously,  N  is  the  number 
of  modes  supported  by  the  acoustic  channel.  For  the  case  of  gaussian  noise,  the 
maximum- likelihood  estimate  of  Vn  for  the  case  of  M  >  N ,  is  given  by  [3]: 

v  =  (M'R-'M)-'M'R-'P,  (5) 

where  the  noise  covariance  matrix  R  is  defined  by  ft  =  E  denotes  expected 

value  and  t  signifies  complex  conjugate  transpose. 

The  present  study  assumes  white  noise  whereby  Eq.  (b)  reduces  to 

\  =  VP  (6) 


with 

V  =  (7) 

This  can  be  seen  by  substituting  ft  =  a2 1  in  Eq.  (5).  Here  /  is  the  identity  matrix. 

Equation  (5)  provides  the  maximum- Ukelihood  solution  for  \,  whereas  we  are  con¬ 
cerned  with  the  estimation  of  r.  The  algebraic  solution  for  r  from  v  however,  is  a 
maximum- Ukelihood  estimator  of  r  if  \  is  the  maximum-Ukelihood  estimate  of  \. 
This  is  guaranteed  by  the  invariance  of  maximum-Ukelihood  estimates  [4], 

A  convenient  form  for  the  solution  and,  in  particular,  for  our  sensitivity  study  can 
be  found  by  defining  a  ‘solution  matrix'  as 


S  =  (vvt),  (8) 

where  the  brackets  indicate  time  average.  Writing  the  pressure  vector  explicitly  in 
terms  of  the  noise-free  signal,  we  have 


P  =  Po  f  f 


(9) 


Thus,  from  Eqs.  (6),  (8)  and  (9)  we  find 

s  =  V(Jw2)r»  4  vrv'.  (io) 
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Equation  (10)  is  derived  based  on  the  assumption  that  the  signal  and  noise  are 
uncorrelated.  Since,  as  pointed  out  earlier,  the  modal  functions  are  real,  the  range 
can  be  computed  from  the  phase  information  alone.  Thus,  from  Eq.  (2c),  one 
obtains 

arg(Sij)  =  (*.  -  *j)r  -  (11) 

where  S, j  is  the  ij th  element  of  the  solution  matrix  S.  The  term  t N  is  necessary 
since  arg( S,;)  only  provides  the  principle  value  of  the  phase  and  the  sign  of  4>„  of 
Eq.  (2a)  is  not  known.  Solving  for  r  yields 


rij 


afg(5u)  ,  y  * 

{ki-h)+  (kt  -  kj) 


(12) 


Defining 


Eq.  (12)  becomes 


0  _  arg(  Si j ) 
r<J  ~  (*,  -  kj)' 


( k-kj )’ 


r'J  =  r?j  +  H&ru. 


(13) 


Since  Eq.  (13)  is  clearly  multiple- valued,  we  must  as  a  minimum  introduce  a  third 
mode.  Upon  introduction  of,  say,  mode  k,  we  would  have  three  solutions:  r,;, 
and  rjk.  Generally,  for  q  modes,  the  number  of  solutions  is  given  by  the  binomial 
coefficient  q,/2,.(q  -  2)!,  i.e.  the  number  of  ways  q  things  can  be  taken  2  at  a  time. 
If  the  solutions  are  computed  out  to  some  maximum  range  and  the  number  of 
solutions  are  counted  up  in  range  bins,  a  plot  can  be  constructed  indicating  the 
number  of  solutions  per  range  bin.  This  is  illustrated  in  Fig.  1.  Here  we  see  the 
number  of  solutions  plotted  for  the  cases  of  4,  6,  8  and  10  modes,  all  of  which 
indicate  25  km  as  the  solution.  This  could  be  thought  of  as  a  pseudo  probability 
density  function.  The  bin  size  here  and  in  all  following  plots  is  100  ni. 


-  3  - 
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3.  Sensitivity  study 


Generally  speaking,  a  sensitivity  study  of  the  type  of  problem  considered  here 
could  ask  two  separate  questions:  (1)  How  sensitive  is  the  algorithm  to  errors  in 
the  measured  field  quantities?  (2)  How  sensitive  is  the  algorithm  to  errors  in  the 
channel  parameters  (model  parameters)?  This  issue  is  discussed  in  detail  ill  the 
appendix.  Here  we  are  interested  in  question  number  (2).  In  particular,  we  will 
look  at  the  sensitivity  to  errors  in  the  SVP  and  the  deviation  of  the  array  from  the 
vertical  (array  tilt).  Also,  we  will  investigate  the  performance  of  the  algorithm  in 
the  presence  of  noise 

The  sensitivity  l<>  errors  in  the  SVP  and  array  angle  were  calculated  based  on  the 
scheme  shown  in  Fig  2.  which  uses  the  SVP  as  an  example.  A  synthetic  data 
sei  is  computed  based  on  the  ’true'  SVP,  i  e.  (SVP)#.  The  signal  that  would  be 
received  on  the  array,  P is  then  computed.  This  synthetic  data  set  is  based  on 
the  modal  functions  designated  as  {<)>„  )0  This  then  is  a  standard  forward  problem. 
Ne«i  a  different  SVP.  designated  (SVP),  is  assumed.  This  differs  from  (S\’  P)o  in 
some  prescribed  wav  and  represents  the  error  introduced  by  our  ignorance  of  the 
true  value  A  new  set  of  modal  functions  is  then  computed.  From  these  modal 
functions,  designated  {</>„},.  the  matrix  1'  (see  Eq  (7))  is  computed.  The  solution 
matrix  ami  the  range  are  then  computed  This  range,  designated  as  r'  in  Fig.  2, 
is  t  h-  n  compared  with  l  the  range  used  in  the  computation  of  the  synthetic  data 
ve,  tor  P„ 


(•or  the  calculations  with  noise,  the  prescribed  noise  matrix  is  used  in  Eq.  (10) 
More  specifically  the  signal  matrix  PP1  was  defined  where 


Ppl 


( s t 

V  N  )  Tr(P,P0')  Tr(RV 


(14) 


where  Tr  designates  trace 

The  scenario  for  the  study  consists  of  an  acoustic  channel  with  a  constant  (in  range) 
depth  of  130  m.  A  point  source  of  190  Hz  is  placed  at  a  depth  of  50  in  and  a  vertical 
array  of  32  equally  spaced  hydrophones  samples  the  acoustic  field  at  ranges  of  10 
and  25  km  This  configuration,  along  with  the  bottom  parameters,  is  summarized 
in  Fig  3  This  corresponds  to  the  area  of  the  Mediterranean  north  of  the  island 
of  Elba  The  bottom  parameters  were  taken  from  (5).  The  sound  velocity  profiles 
were  based  on  measurements  taken  in  the  North  Elba  area  during  the  summer  of 
198fi 

3.1  SENSITIVITY  TO  SVP 

As  mentioned  previously,  the  sound  velocity  data  were  based  on  measurements 
made  during  the  summer  in  the  area  of  interest.  Figure  4  depicts  the  spread  of 
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these  measurements  (slightly  smoothed  for  convenience)  with  the  mean,  which  was 
chosen  as  (SVP)o.  The  parameter  o  is  a  measure  of  the  deviation  of  the  SVP 
(represented  by  the  dotted  curve)  from  (SVP)0.  In  Fig.  4  o  =  0.2  which  means 
that  the  dotted  curve  is  displaced  20%  of  the  distance  between  (SVP)o  and  the 
maximum,  i.e.  the  rightmost  curve.  Thus,  the  dotted  curve  in  Fig.  4  represents 
an  SVP  error  which  constitutes  a  displacement  from  (SVP)o-  This  will  be  referred 
to  as  an  error  of  the  first  type.  A  second  type  of  error  was  investigated.  This  is 
depicted  in  Fig.  5.  Here  the  sign  of  the  deviation  is  alternated,  thus  producing  a 
short-scale  type  error.  This  will  be  referred  to  as  an  error  of  the  second  type.  The 
results  are  shown  in  Figs.  6-12. 

Figures  6  and  7  show  the  results  for  a  range  of  25  km  for  errors  of  the  first  type  for 
a  =  0.1  and  0.2.  These  values  of  a  correspond  roughly  to  sound  speed  deviations 
on  the  order  of  0.5  and  1.0  m/s,  respectively.  Two  things  can  be  seen  from  these 
figures.  First,  the  algorithm  is  quite  sensitive  to  SVP  errors  of  this  type  and  second, 
the  error  impacts  the  higher  order  modes  more  severely.  Turning  to  the  second  type 
of  SVP  error.  Figs.  8  and  9,  which  also  depict  the  25  km  range  case,  indicate  that 
the  results  are  much  less  sensitive  to  errors  of  this  type. 

Here  values  of  a  as  high  as  0.3  (2-3  m/s)  can  be  tolerated  where  for  errors  of  the 
first  type,  values  of  o  greater  than  0.1  cannot  be  tolerated. 

What  is  more  interesting  is  that  errors  of  the  second  type  appear  to  affect  the 
lower-order  modes  more  strongly  than  the  higher-order  modes— quite  the  opposite 
of  the  case  for  the  errors  of  the  first  type. 

Figures  10  and  11  depict  the  results  for  a  range  of  10  km  for  errors  of  the  first 
type  for  a  -  0  1  and  0.5.  Here  we  come  to  the  not  too  surprising  conclusion  that 
the  errors  compound  with  increasing  range.  Also,  it  can  be  seen  that  the  result  is 
biased,  i.e.  the  range  tends  to  be  underestimated. 

Figure  12  depicts  the  10  km  case  for  an  error  of  the  second  type  with  a  -  0.5. 
Again,  we  see  that  the  algorithm  is  less  sensitive  to  this  type  of  error.  As  before, 
we  find  the  conclusion  that  shorter  ranges  can  support  greater  SVP  errors.  Finally 
we  note  that,  as  with  the  25  km  results,  errors  of  the  first  type  tend  mainly  to 
impact  the  higher-order  modes  whereas  errors  of  the  second  type  tend  to  impact 
the  lower-order  modes. 

3.2.  SENSITIVITY  TO  ARRAY  TILT 


Turning  to  the  sensitivity  of  the  algorithm  to  the  tilting  of  the  array  we  refer  first 
to  Figs.  13  and  14.  Here  the  cases  for  ranges  of  10  km  and  25  km  is  depicted 
respectively.  As  can  be  seen,  for  a  tilt  of  0.5°,  the  solution  is  seriously  degraded. 
Beyond  a  tilt  of  1.0°,  the  solution  is  lost. 


-  5 
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The  sensitivity  to  array  tilt  can  be  easily  understood  by  the  following  argument. 
We  begin  by  writing  Eq.  (6)  out  in  explicit  matrix  form: 


'  \l  ' 

Fl.l  El  ,2  1*1 ,32 

Pi 

X.2 

= 

1*2, 1  1*2,2  ■  ■  1*2,32 

Pi 

.  Aio  - 

-  Vio.l  1*10,2  •••  1*10,32  - 

■  P32. 

Here,  we  see  that  each  row  of  the  matrix  V  functions  as  a  beamformer  that  functions 
as  a  spatial  modal  filter. 

In  the  normal-mode  model  each  mode  propagates  as  a  plane  wave  (albeit  with  a 
depth-dependent  amplitude).  Thus,  if  we  apply,  say,  mode  1  to  row  1  of  V  this 
should  extract  mode  1  and  reject  all  the  others.  The  plane  wave  beam  pattern  that 
obtains  when  mode  1  is  applied  to  row  1  is  plotted  in  Pig.  15.  Here  we  see  the 
expected  result  of  maximum  response  at  O°(horizontal).  When  mode  2  is  applied 
to  row  1,  however,  the  beam  pattern  shown  in  Fig.  16  results.  We  see  the  expected 
notch  indicating  the  rejection  of  the  mode.  The  narrowness  of  the  notch  indicates 
that  a  slight  tilting  of  the  array  will  allow  energy  from  the  other  modes  to  spill  in, 
thus  explaining  the  sensitivity  of  the  algorithm  to  the  angle  of  the  array. 


4.  Noise  study 

The  effect  of  noise  on  the  algorithm  was  investigated  for  the  case  of  gaussian  noise, 
both  white  and  coloured.  The  results  for  the  white  noise  case  are  shown  in  Figs.  17 
and  18  which  correspond  to  the  cases  for  a  signal-to-noise  ratio  of  -13  dB.  As 
can  be  seen  from  Eq.  (14),  the  signal-to-noise  ratio  is  that  as  seen  by  a  single 
hydrophone,  i.e.  it  is  the  input  signal-to-noise  ratio.  The  algorithm  broke  down  at 
S/N  =  -15  dB. 

The  results  for  the  case  of  coloured  noise  indicated  no  significant  difference  and 
hence  are  not  included  here. 

Although  these  results  indicate  quite  good  performance  in  the  presence  of  noise, 
they  appear  consistent  with  the  fact  that  the  array  gain  is  of  the  order  of  15  dB. 
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5.  Discussion 


There  are  two  interesting  conclusions  that  can  be  drawn  from  the  SVP  study.  The 
first  is  that  the  results  are  strongly  influenced  by  the  spatial  scale  of  the  SVP  error. 
The  second  conclusion  is  that  which  modes  are  affected  by  the  error  is  strongly 
influenced  by  the  spatial  scale  of  the  SVP  error.  Thus,  based  on  these  results,  one 
sees  that  it  is  the  large-scale  errors  that  are  most  inportant.  This  result  is  quite 
important  since  it  implies  that  making  a  search  of  possible  SVPs  in  the  solution  is 
feasible.  Since  only  the  large-scale  properties  of  the  profile  are  important,  the  class 
of  profiles  that  would  have  to  be  scanned  is  reasonably  small,  since  fine  structure 
is  of  lesser  importance.  In  this  study  the  errors  of  the  second  type  involved  scales 
on  the  order  of  one  wavelength. 

This  insensitivity  to  SVP  fine  structure  suggests  a  possible  extention  of  a  technique 
such  as  this  to  the  active  sonar  case.  With  a  single  active  sonar  ping,  one  can 
determine  the  range  to  a  target.  Given  this  range  and  some  knowledge  of  the 
channel  (i.e.  depth,  bottom  conditions  and  incomplete  SVP  information)  one  could, 
via  such  a  search  method,  conceivably  improve  knowledge  of  the  SVP.  This  suggests 
a  type  of  in-situ  environmental  parameter  update  scheme. 

The  possibility  of  solutions  based  on  parameter  searches  as  discussed  above  also 
presents  itself  in  the  case  of  array  tilt.  Since  the  results  indicate  stability  to  array 
tilts  on  the  order  of  1°,  a  search  of  array  angles  would  present  little  difficulty.  The 
question  of  array  curvature,  however,  is  more  problematical.  Preliminary  studies 
indicate  that  the  algorithm  is  quite  sensitive  to  modest  curvatures.  Dealing  with 
this  problem  by  such  search  techniques  certainly  would  be  more  computationally 
intensive. 

As  mentioned  in  Sect.  4  the  performance  in  the  presence  of  noise  is  quite  good. 
This  is  quite  surprising  in  view  of  the  sensitivity  of  the  algorithm  to  phase  errors 
in  general  (SVP  errors  and  array  tilt).  However,  as  pointed  out  previously,  it  is  not 
inconsistent  if  one  considers  that  the  array  gain  may  be  playing  a  role  here. 
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Appendix  A 

Measurements,  models  and  sensitivity  studies 

In  the  case  of  passive  localization,  the  problem  would  proceed  as  follows.  A  propa¬ 
gation  model  is  selected  (a  normal-mode  model  might  suffice  for  the  shallow-water 
low-frequency  case).  Parameters  required  by  the  model  are  then  obtained  by  mea¬ 
surement  or  from  archival  data.  Measurements  of  the  acoustic  field  are  then  taken 
and  the  position  of  the  acoustic  source  is  estimated  via  some  algorithm.  This  is 
depicted  in  Fig.  Al. 


Channel 


parameters 


Acoustic  field 
measurements 


Source  position 
estimation 


Fig.  At:  Schematic  diagram  of  passive  source  localization  procedure. 


There  are  now  three  questions  that  must  be  asked:  (1)  How  faithful  is  the  model? 
(2)  How  well  must  the  channel  parameters  be  known?  (3)  How  well  must  the  field 
quantities  be  known?  Question  (2),  furthermore,  has  two  parts:  How  accurate  must 
the  measurement  be  and  how  much  time  instability  (fluctuation)  can  be  tolerated? 

Let  Z  be  an  estimated  position  descriptor  (e  g-  range).  Then 

Z  —  Zo  +  AZj i/  +  A  Zp  +  A  Z/,  (dl) 

where  Zo  is  the  true  value,  A  Z\i  is  the  error  due  to  field  measurement  errors,  A  Zp 
is  the  error  in  the  model  parameters,  and  A 2/  is  the  error  due  to  model  infidelity. 

The  accuracy  of  the  estimation  of  Zo  depends  upon  all  three  error  terms  whereas 
the  precision  is  determined  only  by  A  and  A  Zp.  Thus  a  sensitivity  study 
carried  out  with  a  given  model  can  only  determine  the  sensitivity  of  Z  to  errors 
in  the  parameters  and  field  quantities.  That  is,  the  sensitivity  study  can  only 
answer  questions  concerning  precision  (questions  (2)  and  (3))  and  not  accuracy. 
The  accuracy  question  must  be  settled  as  the  second  step  in  the  process.  That  is, 
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given  a  model  a  sensitivity  study  is  carried  out.  If  the  results  indicate  that  the 
problem  is  tractable  (i.e.  A ZM  and  AZp  can  be  rendered  small  enough),  a  series  of 
measurements  must  then  be  performed  to  determine  whether  A Z/  is  small  enough 
to  allow  a  useful  (sufficiently  accurate)  estimation  of  the  range. 

A.l.  SENSITIVITY 

The  sensitivity  problem  can  be  stated  as  follows.  Given  small  changes  in  the  measu¬ 
red  inputs  to  the  model,  what  is  the  corresponding  change  produced  in  the  desired 
quantity?  Thus 


AZM  =  ^  aidti,  ( A2) 

i 

AZp  =  (A3) 

i 

where 


and  {r,  }  and  fy, }  are  the  field  quantities  and  the  model  parameters,  respectively. 
The  sensitivity  problem  then,  is  the  determination  of  the  {a,}  and  {/?,}. 

As  stated  so  far  the  issue  seems  to  be  quite  straightforward.  Via  model  studies,  one 
determines  the  sensitivity  parameters  {a,  }  and  {/?,  }  and  compares  them.  If  the  {«,} 
in  conjunction  with  the  existing  margin  of  uncertainty  in  the  field  measurements 
are  the  major  source  of  error  then  better  measurements  are  called  for.  If,  on  the 
other  hand,  the  {/!,}  indicate  that  the  sensitivity  to  the  model  parameters  is  the 
major  source  of  error,  the  problem  becomes  more  complicated,  since  the  issue  of 
fluctuations  arises.  It  may  simply  not  be  sufficient  to  improve  the  accuracy  of 
the  parameter  measurements  since  the  parameters  themselves  may  be  fluctuating 
in  time.  (The  question  of  spatial  fluctuations  we  consider  to  be  included  in  the 
issue  of  model  fidelity.)  If  this  is  the  case,  and  the  errors  due  to  these  fluctuations 
are  intolerable,  the  only  possible  course  would  be  to  modify  the  model  to  include 
them,  given  that  the  physics  is  known  and  tractable.  It  should  be  noted  here  that 
such  a  course  would  only  guarantee  that  the  model  is  sufficiently  sensitive  and  not 
necessarily  faithful.  That  is,  the  model  may  be  precise  enough  but  still  not  accurate 
enough.  Clearly  though,  this  would  constitute  a  necessary  condition  for  a  faithful 
model.  Since  the  new  model  will  undoubtedly  require  a  new  sensitivity  study,  it 
most  likely  will  follow  that  more  stringent  measurement  requirements  will  arise  for 
the  field  measurements.  Of  course,  this  would  have  been  obvious  from  the  fact 
that  the  fluctuations  in  the  parameters  would  cause  observable  fluctuations  in  the 
measured  field. 

At  first  glance,  the  approach  discussed  above  might  seem  unnecessarily  cumber¬ 
some.  That  is,  why  not  simply  do  the  sensitivity  analysis  in  the  ‘forward  sense’? 
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This  could  be  done  by  varying  the  target  location  in  a  systematic  manner  and  ob¬ 
serving  the  variations  in  the  computed  field.  The  problem  with  this  is  that  it  is 
not  clear  how  such  field  variations  are  related  to  the  efficacy  of  the  algorithm.  This 
problem  is  avoided  in  the  approach  under  discussion  here  since  the  complexities  of 
the  model  and  the  algorithm  are  all  contained  in  the  {a,  }  and  {/?,  }.  It  simply  does 
not  follow  that  if  a  given  error,  say  A Z,  in  the  source  location  produces  a  given 
change  in  phase  or  amplitude  of  the  field  such  phase  or  amplitude  fluctuations  in 
the  measured  field  limit  the  determination  of  the  source  location  to  A Z. 

A  trivial  example  of  this  fact  is  provided  by  the  case  of  a  random  variation  in  the 
amplitude  of  a  narrow-band  source,  where  the  fluctuations  introduced  are  assumed 
to  be  made  up  of  frequencies  much  lower  than  that  of  the  source.  Clearly  in  this 
case,  the  wavenumber  properties  of  the  field,  which  are  the  properties  pertinent  to 
the  spatial  processing  aspects  of  our  problem,  are  essentially  unchanged. 
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Fig.  1:  Plot  of  solutions  for  the  case  of 
25  km  range.  The  number  of  solutions 
per  range  bin  is  plotted  as  a  function 
of  range.  The  bin  size  is  100  m. 


Fig.  2:  Schematic  diagram  of  sensitivity  study  computation. 
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fig  3:  Source-receiver  configuration  and  acoustic  channel  parameters. 
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Fig.  4:  Experimental  SVP  for  the 
region  north  of  Elba.  The  central 
profile  is  the  average  of  a  two- week 
series  of  measurements.  The  outer 
curves  depict  the  spread  of  the 
measurements  and  the  dashed  curve 
represents  a  displacement  of  the 
average  SVP.  In  the  text  this  is  referred 
to  as  an  error  of  the  first  type. 
Here,  the  displacement  is  20%  of  the 
distance  between  the  average  and  the 
maximum,  and  is  denoted  by  a  value  of 
0.2  for  the  parameter  o.  The  numerical 
values  are  those  of  the  mean. 


Sound  Spood  (m) 


1500  1510  1520  1550  1540  1550 


Fig  .*>  SVP  error  of  the  second  type. 
This  curv.  is  identical  to  that  of  Fig.  4 
except  that  the  displacements  alternate 
in  sign. 
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Fig.  8:  Results  for  an  SVP  error  of 
the  second  type:  a  =  0. 1  and  the  true 
range  is  25  km. 
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Fig  9:  Results  for  an  SVP  error  of 
the  second  type:  a  =  0  3  and  the  true 
range  is  25  km. 
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Fig.  10:  Results  for  an  SVP  error  of  the 
first  type:  a  0.1  and  the  true  range 
is  10  km. 


Fig.  11:  Results  for  an  SVP  error  of  the 
first  type:  o  G'5  and  the  true  range 
is  10  km. 
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Fig.  I.’.  Results  for  an  SVP  error  of 
the  second  type:  «  =  0.5  and  the  true 
range  is  JO  km. 


TILT=0.5  DEG. _ M0DES=  1 0 


“Opes- to  ; 

Moots-  a  | 

“■  /  V  ■  A/  /mAv  A.  -  ./A 

/Av4^-  \ma 

-  soots-.  6  1 

•saJV/VV'^A >/y\  /  vA  >A/  /  vAv 

MOOtS-  *  J 

Range  (km) 


Fig  13:  Results  for  an  array  tilt  of  0.5”: 
the  true  range  is  10  km. 
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Fig  14:  Results  for  an  array  tilt  of  0.5°: 
the  true  range  is  25  km. 


Fig.  15:  Angular  response  of  the 
vertical  array  to  mode  no.  1:  the 
beamfornrer  weights  are  the  modal 
filter  coefficients  fo’r  mode  no.  1. 
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Fig.  16:  Angular  response  of  the 
vertical  array  to  mode  no.  2:  the 
beaniformer  weights  are  the  modal 
filter  coefficients  for  mode  no.  1. 
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Fig.  17  Results  for  the  case  of  white 
gaussiau  noise:  the  signal-to- noise 
ratio  is  -13  dB  and  the  true  range  is 
10  km 
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